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A method of determining the sound pressure resulting from a surface 
element of a sound emitting surface 

Fieid of me invention 

The Invention relates to the detection and identification of sound sources, 
and noise sources in particular, with the purpose of reducing noise emission. 

Background of the invention 

Protection of the environment and human beings has become increasingly 
important. Buildings, vehicles such as cars, buses and aircraft, household 
appliances and industrial machinery have noise producing components such 
as engines, motors, gears, transmissions etc. In order to protect individuals 
from such noise, the noise generating components and the transmission path 
of the noise to a human being have been investigated with the purpose of 
reducing the generated noise at the source and of reducing the noise trans- 
mitted from the source to human beings. 

Testing of acoustic properties of noise generating and noise transmitting me- 
dia such as mechanical structures and air or other fluids is an important part 
of the process of noise reduction. In complex structures with several noise 
sources, such as mentioned above, it can be complicated to identify noise 
sources and transmission paths and their contributions to the perceived 
noise. 

Mathematical models and computerised methods exist for vibro-acoustic 
analysis of physical structures. Acoustical tools exist for simulating acoustic 
properties of portions of a human being, such as Mouth Simulator type 4227, 
Ear Simulators types 4185 and 4195, Head and Torso Simulators types 4100 
and 4128. all from Bruel & Kjaer Sound and Vibration Measurement A/S, 
Denmark. All of these are intended for use in analysing sound at different 
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stages in its "normal" forward transmission from the source to a human be- 
ing. 

EP 0 015 852 discloses a three-dimensional array of microphones for meas- 
uring the total or directional acoustic power emitted by a sound source. Such 
5 an array is suitable for use with the method of the present invention. 

The transfer function for sound from a small omni-directional sound source to 
a point of measurement is often expressed as the acoustic transfer function 
H (or acoustic transfer impedance Z t ) defined as H = p/Q, where Q is the vol- 
ume velocity emitted from the sound source, and p is the sound pressure at 

10 the point of measurement resulting from the volume velocity Q generated by 
the sound source. In most cases the analysed mechanical and acoustic 
transmission media are reciprocal, which means that the acoustic transfer 
function is the same both for forward and reverse transmission. In other 
words, if the sound source and the measuring microphone are interchanged, 

15 whereby the transmission of sound through the media is reversed, and 
boundary conditions remain unchanged, then the acoustic transfer imped- 
ance is unaffected, ie the "forward" acoustic transfer impedance and the "re- 
verse" acoustic transfer impedance are identical. 

It is known to use this fact when analysing the transmission of sound, 
20 whereby a sound source is placed in a position that is normally occupied by a 
human being, ie a "listening" position, and a microphone is placed in the 
normal position of the sound source. This has distinct advantages when iden- 
tifying sound sources and tracking the noise on its path from the source to 
the listening position. 

25 When measuring the forward transmission path a Head and Torso Simulator 
type 4100 from Bruel & Kjaer Sound and Vibration Measurement A/S can be 
placed in the listening position, whereby very realistic measurements of the 
forward transmission path can be obtained, since the influence of the head 
and the torso on the transfer function to the ears is taken into account Dan- 
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ish patent application PA 200300589 discloses a simulator simulating acous- 
tic properties of the head and possibly the torso of a human being. That 
simulator comprises a sound source for outputting sound signals through the 
simulated ears. Such a simulator completes the reverse measuring chain and 
5 can be placed in a position that is normally occupied by a human being, ie a 
listening position 1 *. By means of a pair of microphones in each simulated ear 
canal the output sound volume velocity can be measured. This is useful for 
computing the acoustical transfer function from a sound source to a listening 
position. 

10 When designing eg vehicles such as cars, buses and aircraft the comfort of 
the passengers, the driver and crewmembers is of importance. Noise can 
seriously jeopardize not only comfort but also the health of humans. It is 
therefore important to reduce noise ( and for effectively reducing noise it is 
important to identify noise sources and their individual contribution to the 

15 noise level at locations where people are present. Mechanical structures 
such as body and wall panels can vibrate and emit noise, and large struc- 
tures can have "hot spots" that emit more noise than "cold spots". Not all hot 
spots may be serious contributors to the noise level resulting at a "listening 
position", and, vice versa, cold spots may contribute more seriously than ex- 

20 pected. Such phenomena can e.g. be due to conditions in the transmission 
path from the source to the listening position. 

The problem to be solved by the invention is to provide a method of deter- 
mining, in a predefined position such as a listening position of a human be- 
ing, the sound pressure resulting from sound emitted from a surface element 
25 of a sound emitting surface. In particular there is a need for identifying, 
among the plurality of noise sources, which can be distributed over a large 
area, the most significant sources and their contributions to the noise level at 
one or more listening positions. 
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Summary of the invention 

The invention provides such a method for determining the contribution of 
each surface element of a sound-emitting surface to the noise level at e.g. a 
listening position. 

5 The sound pressure in the listening position resulting from the sound emitted 
from a surface element can be calculated by multiplying the emitted volume 
velocity by the acoustic transfer function (also referred to as the acoustic 
transfer impedance) from the surface element to the listening position. In 
general, the acoustic transfer function is defined as the (complex) ratio of the 
10 effective sound pressure p at a given point to the volume velocity q of the 
sound source generating the sound pressure p. 

It is assumed that the analysed mechanical and acoustic transmission media 
are reciprocal, which means that the acoustic transfer function is the same 
both for forward and reverse transmission. The acoustic transfer function can 
1 5 advantageously be found as the "reverse" acoustic transfer function. 

The invention comprises the following two main phases, in which two differ- 
ent measurements are performed to determine the noise source strength and 
the transfer function, respectively. 

I. An "operational measurement" is made with the noise source(s) operating 
20 under the conditions to be investigated. A three-dimensional array of micro- 
phones is used to measure the three-dimensional sound field above a sur- 
face element of a sound emitting surface, and at the surface element the 
component of the air-particle velocity perpendicular to the surface element is 
calculated using e.g. Near-Field Acoustical Holography (NAH), the Statisti- 
25 cally Optimal NAH (SONAH) as described in [2], or the Inverse BEM (Bound- 
ary Element Modelling). The volume velocity emitted by the surface element 
is then determined by integrating the air-particle velocity over the area of the 
element - or approximately as the air-particle velocity perpendicular to the 
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surface element multiplied by the area of the surface element. Several sur- 
face elements may be covered and measured simultaneously, and the di- 
mensions of each element should be small in comparison to the wavelength. 

II. A "transfer function measurement" is made with the noise source(s) under 
5 investigation inactivated. This measurement is typically taken with the same 
array of microphones in the same position as in the operational measurement 
above, but now the surface element itself must be rigid and non-vibrating. A 
volume velocity sound source is used to emit a volume velocity in a listening 
position, and the array of microphones is used to determine the resulting 

10 three-dimensional sound field above the surface element. The sound emitted 
from the volume velocity sound source should preferably be the only sound 
or at least dominate over possible other sounds. Using e.g. Near-Field 
Acoustical Holography the resulting sound pressure at the surface element is 
determined. The acoustic transfer function between the surface element and 

15 the listening position is assumed to be reciprocal and is determined as the 
ratio of the resulting sound pressure at the surface element to the volume 
velocity emitted by the volume velocity sound source. 

The sound pressure in the listening position resulting from sound emitted 
from the surface element is then determined as the volume velocity emitted 
20 by the surface element multiplied by the acoustic transfer function. 

For each position of the microphone array 10 all the above steps are re- 
peated, and a noise source "map" can be made of the entire interior surface 
of e.g. a car cabin, or important parts thereof. Such a map can be used for 
identifying the most serious noise sources to be attenuated. 

25 The above-described method requires the measurement of the reciprocal 
transfer function to be performed with the air-particle velocity perpendicular to 
the surface elements equal to zero, which in turn requires that the surface 
elements are rigid. According to the above, the only acoustical quantities on 
the surface element that are taken into account are the "operational" air- 
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particle velocity at the surface and the sound pressure resulting from the vol- 
ume velocity sound source. 

A more general method uses also the "operational 11 sound pressure and the 
particle velocity created by the volume velocity sound source, whereby the 
5 requirement for rigid surface elements is avoided. 

Brief description of the drawings 

Figure 1 is a plan view of a portion of a three-dimensional array of micro- 
phones used with the invention, 

Figure 2 is a side view of the three-dimensional array of microphones in fig- 
10 ure 1 shown in a measurement relationship to a sound-emitting surface, and 

Figure 3 illustrates a measurement set-up for measuring in e.g. a car cabin. 

Detailed description of the invention 

The measurement set-up used in the invention is as follows. 

in figures 1 and 2 a three-dimensional array 10 of a plurality of microphones 
15 M is shown. The term "microphone" is here used as the generally accepted 
term for a transducer that generates an (electrical) output signal in response 
to pressure variations in the fluid (air, water etc.) that surrounds the trans- 
ducer The array comprises two identical, parallel layers of microphones M 
spaced a vertical distance d v apart In each layer the microphones are dis- 
20 tributed in two sets of parallel rows forming a square grid with the rows 
spaced horizontal distances dhx and dhy, respectively, apart. In the preferred 
embodiment shown, the vertical distance d v and the horizontal distances dhx 
and dhy are identical, whereby the microphones are uniformly distributed and 
form a cubic lattice. A non-uniform distribution of the microphones can be 
25 used. The microphones have well-defined, preferably identical, electro- 
acoustical properties. The array of microphones preferably has a multi-pole 
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plug (not shown) for connecting the microphones to measuring equipment, 
which is not part of the invention. Each layer of the array can have eg 6x6 or 
8x8 or any other suitable arrangement of microphones. The vertical and hori- 
zontal spacing determine the upper frequency limit at which the array can be 
5 used. A vertical and horizontal spacing of 5 cm results in an upper frequency 
limit of about 3 kHz. 

As indicated in figure 2, the two layers are separate layers that are assem- 
bled and can be disassembled and used independent of each other. The two 
layers are here planar layers mounted back-to-back, but a fixed three- 
10 dimensional array may also be used. 

In figure 3 is shown a rectangular frame illustrating e.g. the cabin of a car. 
The cabin has an interior surface S, which, due to vibrations from the engine, 
the tyres etc that are transmitted to the cabin, will emit corresponding sound 
(noise) into the cabin. 

15 In figures 2 and 3 the three-dimensional array 10 of microphones is arranged 
above the surface S. A surface element AS of the surface S is shown be- 
neath the microphone array 10. The position of the microphone array 10 rela- 
tive to the surface element AS including the distance D and the lateral coor- 
dinates must be well-defined. For regular array geometry, there will typically 

20 be one surface element AS for each microphone (in the layer close to the 
surface), with the microphone at the centre of the corresponding surface 
element. 

With the above set-up the following measurements are taken. 

First, an "operationar measurement Is taken with the set-up shown in figures 
25 2 and 3. With the "real life" noise sources (engine, tyre noise etc) active, and 
with the microphone array 10 in the shown position, the sound pressure at 
each of the microphones M is measured. 
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Based on the measured three-dimensional distribution of sound pressure 
over the array of microphones, a calculation is performed of the air-particle 
velocity u n perpendicular to the surface element AS resulting from the sound 
emitted from the entire surface S. This can be done using e.g. a well-known 
5 Near-Field Acoustical Holography (NAH) method, whereby the three- 
dimensional acoustical near-field is extrapolated to the surface element AS. 

Next, a "transfer function measurement" is made. A volume velocity sound 
source is used to emit a well-defined volume velocity at a predefined position 
such as a listening position, i.e. a position In which a person occupying a seat 
in the car would normally have his head or one of his ears. A volume velocity 
sound source is a sound source that has a speaker in an enclosure with an 
orifice through which the speaker emits sound with a well-defined volume 
velocity. The output volume velocity can be measured/calculated by well- 
known methods, or the volume velocity sound source can be calibrated to 
emit a volume velocity that is known due to calibration. 

The preferred volume velocity sound source 1 1 is the one disclosed in Dan- 
ish patent application PA 200300589, which is a simulator simulating acous- 
tic properties of the head and the torso of a human being. The simulator has 
a sound source for outputting well-defined volume velocity sound signals 
20 through the simulated ears for reciprocal measurement of the transfer im- 
pedance from a source position to one of the simulated ears. 

The noise sources, the engine etc., are switched off, and the volume velocity 
sound source 11 is used as at least the dominating sound source and pref- 
erably the only sound source to avoid possible influence from other sound 
25 sources. With the three-dimensional array of microphones still in a predefined 
position (preferably the same as above) relative to the surface element, the 
sound pressure at each of the microphones, Le. the three-dimensional distri- 
bution of sound pressure over the array of microphones, is again measured. 
Using e.g. a Near-Field Acoustical Holography (NAH) method, the three- 
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dimensional distribution of sound pressure over the array of microphones, i.e. 
the acoustical near-field is extrapolated to calculate the sound pressure on 
the surface element AS resulting from the sound emitted by the volume ve- 
locity sound source. 

5 The transfer function is then calculated as the ratio of the calculated sound 
pressure at the surface element to the volume velocity emitted from the vol- 
ume velocity sound source. The thus calculated transfer function is both the 
reverse transfer function and the forward transfer function. 

With the volume velocity emitted by the surface element and the transfer 
10 function both being known, the resulting sound pressure at the listening posi- 
tion can be found by multiplying the volume velocity by the transfer function. 

Depending e.g. on the size of the microphone array 10 and the size of the 
surface element AS, measurement and calculation of the sound pressure 
resulting from one or more surface elements AS can be performed with one 
15 position of the microphone array 10. If the contributions from several surface 
elements are calculated, then these can be added to obtain the contributions 
from larger areas. 

The above-described method requires that the air-particle velocity u„ perpen- 
dicular to the surface element AS is zero when measuring the transfer func- 
20 tion. This requirement is only fully satisfied when the surface S is rigid. In 
case of a car cabin the interior surfaces often have sound attenuating or 
sound absorbing materials, which are not rigid, or such materials delimit the 
useful space of the cabin. 

The more general version of the method of the invention to be described be- 
25 low is not limited to rigid surfaces, but the surfaces can be soft such as sound 
attenuating or absorbing materials. On soft materials the air-particle velocity 
u n perpendicular to the surface element AS is not necessarily zero when 
measuring the transfer function, i.e. with sound emitted from a volume veloc- 
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ity sound source, fn the more general version of the method of the invention 
the same measurements as described above are taken, but the processing of 
the measurements is slightly extended with more calculations. 

Like above described, an "operational" measurement is taken with the set-up 
5 shown in figures 2 and 3. With the noise sources (engine, tyre noise etc) ac- 
tive, and with the microphone array 10 in the shown position, the sound pres- 
sure at each of the microphones M is measured. In addition to the calculation 
of the air-particle velocity u n perpendicular to the surface element AS result- 
ing from the sound emitted from the surface element, the sound pressure p 
10 on the surface element AS is also calculated. This, too, can be done using 
e.g. a well-known Near-Field Acoustical Holography (NAH) method. 

Further, also a "transfer function measurement" is made. This is done with 
the set-up shown in figure 3. In addition to the calculation of the sound pres- 
sure pv on the surface element AS resulting from the sound emitted by the 
15 volume velocity sound source, the air-particle velocity u v .n perpendicular to 
the surface element AS resulting from the sound emitted by the volume ve- 
locity sound source is also calculated. This, too, can be done using e.g. a 
well-known Near-Field Acoustical Holography (NAH) method. 

The sound pressure Ap at the listening position is then calculated using the 
20 following formula 

where Q v is the volume velocity emitted by the volume velocity sound source. 

For practical purposes, when the variations in sound pressure and particle 
velocity over the surface element AS are small and can be regarded as con- 
25 stant, the surface integral in the above formula (1) can be calculated using 
the following more simple formula 
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The method of the invention is not restricted to dosed volumes such as 3 car 
cabin, but can also be used for determining the sound pressure resulting 
from virtually any distributed sound source such as large machinery or a plu- 
5 rality of machines. 

The surface element AS needs not be the vibrating surface of a mechanical 
device but can be any imaginary surface in an acoustic medium such as air, 
or the surface element AS may include or be part of a sound-emitting open- 
ing in a mechanical device. 
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Claims 

1 . A method of determining, in a predefined target position, the sound pres- 
sure (Ap) resulting from sound emitted from a surface element (AS) of a 
sound emitting surface (S), the method comprising 

5 - measuring, using a three-dimensional array of a plurality of microphones 
arranged in a first predefined measuring position relative to the surface ele- 
ment (AS), a first three-dimensional distribution of sound pressure, and 

- calculating, based on the first three-dimensional distribution of sound pres- 
sure, the air-particle velocity (u n ) on the surface element (AS) and perpen- 

10 dicular to the surface element (AS), resulting from the sound emitted from the 
surface (S), 

- arranging a sound source capable of emitting a volume velocity (Q v ) in the 
target position, 

- causing the sound source to emit the volume velocity (Qv). 

15 - measuring, using a three-dimensional array of a plurality of microphones 
arranged in a second predefined measuring position relative to the surface 
element (AS) and with the volume velocity (Q v ) emitted from the sound 
source in the target position creating a dominating sound, a second three- 
dimensional distribution of sound pressure, and 

20 - calculating, based on the second three-dimensional distribution of sound 
pressure, the sound pressure (pv) at the surface element (AS) resulting from 
the volume velocity (Q v ) emitted from the sound source in the target position, 

- determining the transfer function H = p«/Q v as the ratio of the sound pres- 
sure (p v ) at the surface element (AS) to the volume velocity (Q v ) emitted from 

25 the sound source in the target position, and 
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- determining the sound pressure (Ap) in the target position as 

AD = H./lu-ASV 

2. A method of determining, in a predefined target position, the sound pres- 
sure (Ap) resulting from sound emitted from a surface element (AS) of a 
5 sound emitting surface (S), the method comprising 

- measuring, using a three-dimensional array of a plurality of microphones 
arranged in a first predefined measuring position relative to the surface ele- 
ment (AS), a first three-dimensional distribution of sound pressure, and 

- calculating, based on the first three-dimensional distribution of sound pres- 
10 sure, the air-particle velocity (u n ) perpendicular to the surface element (AS) 

and on the surface element (AS), and the sound pressure (p) on the surface 
element (AS), resulting from the sound emitted from the surface (S), 

- arranging a sound source capable of emitting a volume velocity (Q v ) in the 
target position, 

15 - causing the sound source to emit the volume velocity (Q v ), 

- measuring, using a three-dimensional array of a plurality of microphones 
arranged in a second predefined measuring position relative to the surface 
element (AS) and with the volume velocity (Q v ) emitted from the sound 
source in the target position creating a dominating sound, a second three- 

20 dimensional distribution of sound pressure, and 

- calculating, based on the second three-dimensional distribution of sound 
pressure, the sound pressure (pv) at the surface element (AS) and the com- 
ponent of the particle velocity (u v , n ) perpendicular to the surface element (AS) 
resulting from the volume velocity (Q v ) emitted from the sound source in the 

25 target position, 
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- determining the sound pressure (Ap) in the target position in accordance 
with the formula 
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3. A method according to any one of claims 1-2 characterized 
5 in that the target position is a listening position suitable for being occupied by 

a human being. 

4. A method according to claim 1 characterized in that the 
air-particle velocity (u n ) perpendicular to the surface element (AS) resulting 
from the sound emitted from the surface (S) is calculated, based on the first 

10 three-dimensional distribution of sound pressure, using a Near-Field Acousti- 
cal Holography (NAH) method, and that 

the sound pressure (pv) at the surface element (AS) resulting from the vol- 
ume velocity (Q v ) emitted from the sound source in the target position is cal- 
culated, based on the second three-dimensional distribution of sound pres- 
15 sure, using a Near-Field Acoustical Holography (NAH) method. 

5. A method according to claim 2 characterized in that the 
air-particle velocity (u n ) perpendicular to the surface element (AS) and the 
sound pressure (p) resulting from the sound emitted from the surface (S) are 
calculated, based on the first three-dimensional distribution of sound pres- 

20 sure, using a Near-Field Acoustical Holography (NAH) method, and that 

the sound pressure (p v ) at the surface element (AS) and the air-particle ve- 
locity (uv, n ) perpendicular to the surface element AS resulting from the vol- 
ume velocity (Q v ) emitted from the sound source in the target position are 
calculated, based on the second three-dimensional distribution of sound 
25 pressure, using a Near-Field Acoustical Holography (NAH) method. 
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6. A method according to any one of claims 1-5 characterized 
by using as the volume velocity sound source a simulator simulating acoustic 
properties of at least a head of a human being, the simulator having a simu- 
lated ear with an orifice and a so unci source for uuiputting sound slgnais 

5 through the orifice of the simulated ear. 

7. A method according to claim 6 characterized in that the 
simulator simulates the acoustic properties of the head and a torso of a hu- 
man being. 

8. A method according to any one of claims 1-? characterized 
10 by using, as the three-dimensional array of a plurality of microphones, an 

array having two parallel layers of microphones, where each layer comprises 
a plurality of microphones arranged in a two-dimensional grid. 

9. A method according to claim 2 characterized in that the 
sound pressure (Ap) in the target position is determined as an approximation 

1 5 in accordance with the formula 
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A method of determining the sound pressure resulting from a surface 
element of a sound emitting surface 

Abstract 

A three-dimensional array (10) of microphones (M) is used to determine the 
5 three-dimensional sound field above a surface element (AS) of a sound emit- 
ting surface (S), and the air-particle velocity (u n ) at the surface element (AS) 
is determined using Near-Field Acoustical Holography (NAH). A volume ve- 
locity sound source (11) is used to emit a volume velocity (Q v ) in a listening 
position, and the array (10) of microphones (M) is used to determine the re- 

10 suiting three-dimensional sound field above the surface element (AS), and 
using NAH the resulting sound pressure at the surface element (AS) is de- 
termined. The acoustic transfer function (H) between the surface element 
(AS) and the listening position is assumed to be reciprocal and is determined 
as the ratio of the resulting sound pressure at the surface element (AS) to the 

15 volume velocity (Q v ). The sound pressure in the listening position resulting 
from the surface element (AS) is determined as Ap = H(u n AS). 

Figure 3 should be published 
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